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ABSTRACT 

Motivated by the discovery of a number of radio relics we investigate the fate of fossil 
radio plasma during a merger of clusters of galaxies using cosmological smoothed- 
particle hydrodynamics simulations. Radio relics are extended, steep-spectrum radio 
sources that do not seem to be associated with a host galaxy. One proposed scenario 
whereby these relics form is through the compression of fossil radio plasma during 
a merger between clusters. The ensuing compression of the plasma can lead to a 
substantial increase in synchrotron luminosity and this appears as a radio relic. Our 
simulations show that relics are most likely to be found at the periphery of the cluster 
at the positions of the outgoing merger shock waves. Relics are expected to be very 
rare in the centre of the cluster where the life time of relativistic electrons is short and 
shock waves are weaker than in the cooler, peripheral regions of the cluster. These 
predictions can soon be tested with upcoming low-frequency radio telescopes. 

Key words: galaxies: clusters: general - intergalactic medium - shock waves, cos- 
mology: theory - diffuse radiation, radiation mechanism: non-thermal 



1 INTRODUCTION 

Active galactic nuclei (AGN) inject a large amount of mag- 
netised, relativistic plasma into the intra-cluster medium 
(ICM). This radio plasma emits mainly synchrotron ra- 
diation. However, after a typical time of 10* years the 
plasma has cooled radiatively such that the rem aining radio 
emission is difficult to detect fe.g. Ijaf^ I^T^)). The rem- 
nants of radio lobes are called 'radio ghosts' or 'fossil radio 
plasma'. Recently, substantial evidence for radio ghosts has 
come from the detection of cavities in X-ray surface bright- 
ness maps of clusters of galaxies fBoehringer et al]|l99j; 
fchurazov et al. 200 0; Wils on ct al. 2000: McNamara et alJ 
BoOlinBlanton et alJr200lh . 



Diffuse, steep-spectrum radio sources with no optical 
identification have been observed in a growing number of 
galaxy clusters. These objects have complex morphologies 
that show diffuse and irregular emission in combination 
with point-like sources. They are usually subdivided into 
two classes: those that are located near the centre of a 
clust er, e.g. in Abell 520 and Abell 2254 (Giovannini ct al. 
I1999I) . and those that are located in the p eriphery of a clus- 
ter, e . g. Abell 85 Abe ll 133, Abell 3667 llRottgering et alJ 
I1997I: Isiee et aljEool . denoted as 'radio halos' and 'ra- 
dio relics', respectively. Unlike halos, radio relics have a 
filamentary morphology and show a partial polarisation 
of the radio emission. This distinction, however, is not 
free of contradictions. The cluster Abell 520, e.g., shows 
knotty radio structures located in the centre of the clus- 
ter's X-ray emission. The rough classification reflects the 



poor understanding of the origin of the radio sources and 
raises the question whether halos and relics are indeed pro- 
duced by different processes. R adio halos and relics show 
a steep spectrum a 1 — 1.8 iKempner fc Sarazi nl l200ll : 
iBacchi et alll2003ll and the cut-off at high freq uencies in- 
dicates that the electron p opulation has aged iSlee et alJ 
l200lHKaiser fc CotteJliooil . For more details on observa- 
tions of diffuse cluster radio sources the read er is referred to 
jKempner fc SarazinllioOll : lGiovanninilll999l) . The difliculty 
in explaining the radio emission lies in the lack of an evident 
source for the relativistic electrons, such as an AGN. The 
strongest hint for the formation of, both, halos and relics 
may come from the fact that both are observed in clusters 
that show signs of an ongoing merger, for instance signifi- 
cant substruc tures in the X -ray emission or the absence of a 
cooling fiow jBacch^^^ HSP'^ • Kcmoncr fc Sarazin 20oJ 
Ferett i 'l999'; Schuecke rfc Bohringed Il999l: IVenturi et all 
19991: iRoettiger et al..>1999r " "n the rest of the paper we will 
be primarily concerned with radio relics because their for- 
mation is believed to be quite different from that of radio 
halos. 

Shock waves that are produced by a merger between 
clusters of galaxies may provide the necessary acceleration 
of the electrons. Several processes for the formation of radio 
relics have been proposed, the two most important being: 
(i) i n-situ diffusive shock acceleration by the Fermi I pro- 
cess jEnsslin et aljl998l : IRoettiger et aljl9'99l : iMiniati et alJ 
I2OOII) and (ii) re-acceleration of el ectrons by compression of 
existing cocoons of radio plasma (lEnfilin fc Gopal-Krishnal 
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Figure 1. The density (upper panel) and temperature (lower panel) distribution on a slice along a plane which contains the both centres 
of mass of the initial two clusters. The length scale is given in comoving Mpc. The first snapshot is taken at redshift z = 0.66 and 
the later snapshots follow approximately the motion of centre of mass of the cluster. The density is given in proper mass density divided 
by the proton mass. 



I2OOII) . In all the aforementioned scenarios the radio emis- 
sion traces the shock-front. The fact that radio relics re- 
semble individual objects and do not trace the entire shock 
front, may provide some indication in favour of the last sce- 
nario. Moreover, when a radio ghost is passed by a cluster 
merger shock wave with a typical velocity of a few 1000 
km/s the ghost is compressed adiabatically and not shocked 
because of the much higher sound speed within it. There- 
fore, diffusive shock acceleration is unlikely to be the prime 
mechanism that re-energises the relativistic electron popu- 
lation. However, it has been shown that the energy gained 
during the adiabatic compression together with the increase 
in the magnetic fields strength c an cause the fossil radio 
cocoon to emit radio waves again. lEnfilin fc Go pal-Kris hnal 
l|200 J) showed that the spectra thus produced are consistent 
with an old electron population that has been adiabatically 
compr essed. With the aid o f mag neto-hydrodynamical simu- 
lations ^nfih^^^Brtgge^ (|2002|) demonstrated that the re- 
sulting relics possess a toroidal shape and are partially polar- 
ized - in good agreement with observations. One prerequisite 
for this mechanism to be effective is that the electron popu- 
lation is not older than 0.2 - 2 Gyr (Enf51in & Gopal-Krishna 
2001). The time-scale depends on the conditions in the sur- 
roundings, mainly the external pressure. In high-pressure en- 
vironments, such as in cluster cores, the synchrotron losses 
are expected to be higher than in regions of lower pressure 
because the magnetic fields in the pressure-confined radio 
plasma is higher. This leads to a shorter synchrotron cool- 
ing time for the fossil radio plasma and a reduced probability 
to fiare up during the passage of a shock. 

In this paper we use numerical simulations to study the 
ICM during a merger of galaxy clusters. Here we focus on the 
formation of radio relics and study, in particular, whether 



shock waves, that are associated with the merger, can re- 
vive the fossil plasma. The main issues that we would like 
to address are (i) if radio cocoons that have been produced 
before the onset of the merger can be revived by shocks, 
and (ii) where the revived radio plasma is most likely found. 
To produce radio relics, we apply the formalism suggested 
bv lEnfjlin fc Gopal-Krishnal ll200 j) to high-resolution simu- 
lations of a major merger. 



2 SIMULATIONS 

We performed high-resolution simulations of mergers be- 
tween clusters of galaxies using the Gadget dSpringel et alJ 
I2OOII) code, which is a smoothed-particle hydrodynamics 
(SPH) code based on a tree-scheme to compute the grav- 
itational interaction. We used the entropy conserving SPH- 
kernel as proposed by ISpringel fc Hernauisd (|2003) who 
found that using the entropy as independent variable results 
in an improved representation of explosion shock waves. 

Gadget provides the possibility to use particles with 
different masses. This allows us to simulate a cluster merger 
with a high mass-re solution following the scheme given by 
iKlvpin et alJ i200ll) . The simulations are set up according 
to the standard ACDM model with 57a/ = 1 — f^A = 0.3, 
Q.B ~ 0.039, and h = 0.7, where the Hubble constant is given 
by 100 /ikm, s~^ Mpc~^. For redshift z — 50 the initial par- 
ticle positions and velocities are derived using a power spec- 
trum normalised on the 8/i~^Mpc scale to erg = 0.9. Initially, 
a simulation in a computational box of 80h^^ Mpc with 128'' 
particles is performed. When the simulation arrives at z — 0, 
a region with a suitable cluster is isolated. In a new simula- 
tion the mass-resolution is refined in this region by splitting 
the particles up into 64 sub-particles. Thus, we achieve a 
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Figure 2. The radio luminosity probability at a frequency of 100 MHz per mass of radio plasma in the same slice through the cluster 
centre as in Fig. The prescription for the evolution of the relativistic electron distribution is described in Sec. |21 Here we assumed a 
pressure ratio between the magnetic and gas pressure of Pb/P^ss = 0.1 and 0.01 in the upper and lower panels, respectively. 



gas mass-resol ution of 3.6 x 10'^ h ^^ Mq in the cluster re- 
gion. Following lPovyer et al.l l)2nn3|) we choose a gravitational 
softening length for the refined region 

nal cluster is produced by a merger between two progenitors 
with similar mass, 1.7x lO" Mq and 1.6 x 10^^ Mq, 
at redshift z ~ 0.6. The collision processes virtually central 
with a relative speed of 2000 kms^^. 



We should note that effects of gas cooling, stellar feed- 
back, thermal conduction, and magnetic fields are not in- 
cluded in this simulation. Our primary goal is to model the 
evolution of the merger shock waves in the ICM within a 
realistic cosmological setting. In order to calculate the radio 
emissivity of non-thermal plasma, we have to make assump- 
tions about the strength of the magnetic field. Very little is 
known on the strength of the magnetic field in radio ghosts. 
Measurements of the large-sclae mag netic field in clusters in - 
dicate a field strength of few /LtGauss iCarilli fc Tavloil2002l) . 
In the centre of the cluster the field strength is higher and ap- 
proximately thermal, i.e. B^/Stt = Pb ~ Pgas, (Eilck 199£|. 
However, it seems unlikely that there is enough small-scale 
turbulence or that there are enough powerful radio galax- 
ies to sustain a thermal magnetic field throughout the en- 
tire cluster volume. Moreover, these measurements rely on 
Faraday rotation measures, to which the dilute plasma in 
the radio ghosts does not contribute significantly. Therefore, 
these measurements can only serve as a very rough guide to 
the actual field strength in the bubbles. In our calculations 
we assume that the fraction between magnetic pressure and 



thermal pressure is fixed at a value which yields a mass- 
averaged field strength of a few /xGauss. 



3 THE FATE OF RADIO HALOS 

Relativistic electrons in radio lobes lose energy via syn- 
chrotron radiation and inverse Compton scattering with 
photons of the cosmic microwave background (CMB). In 
addition, adiabatic compression of the cocoons is able to re- 
accelerate the electrons. On the basis of our simulation we 
investigate whether the compression that the radio plasma 
suffers in the course of the merger can compensate for its 
radiative losses and is thus able to revive the fossil radio 
plasma. In the following, we assume that each gas particle 
also carries a certain amount of non-thermal radio plasma 
with it, whose evolution is calculated as follows. 

The momentum of a relativistic electron in a radio lobe 
is altered by synchrotron losses that are proportional to the 
magnetic energy density iis, by inverse Compton losses pro- 
portional to the CMB field density ucmb and by adiabatic 
compression of the considered volume V 

4 11 
dp = --CTT {ub + ucmb} - -^y dV , (1) 

where ctt denotes t he T homson cross section. 
lEnfihn fc Gopal-Krishiial J200ll) showed that the the 
electron distribution function f{p,t) dp dV can be derived 
from the initial distribution /(p, to), the compression ratio 



Cit) = Vito)/V{t) = (P(t)/P(to))^/^ : 



(2) 
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where the volume V of the radio ghosts is adiabatically com- 
pressed by the ambient pressure P with an adiabatic expo- 
nent 7 for magnetized plasma, and the characteristic mo- 
mentum 

i- = lax dt' {uB{t') + ucMB(t')} (^) • (3) 

If the initial distribution is a power-law f{p,to) = 
fo(p/po)~" the spectrum at later times becomes 

f{p,t) = /oc(t)<'-"'/^ ip/vo)-" {i^p/p*r-\ (4) 

The luminosity of a fossil radio plasma depends thus on the 
age of the radio plasma, the strength of the magnetic field, 
the CMB density and the pressure history of the ambient 
gas. To set a starting time io of the pressure history we 
assume that the radio plasma has been released at an early 
stage of the merger, more precisely when the two progenitors 
are still separated by 0.5 Mpc. Those ghosts may result 
from the advection of the radio lobes by the merger, i.e. 
the merger may itself have generated the radio ghosts. The 
assumption that the radio plasma in the ghosts dates from 
the beginning of the merger yields a conservative estimate 
for the later occurrence of radio relics. As argued above, 
the magnetic field in the outer part of the cluster should 
be sub-thermal. We assume that the ratio Pb / Pgas between 
the magnetic and thermal pressure is in the range of 1-10%. 
Furthermore, we assume that the fossil radio plasma moves 
with gas, i.e. the plasma does not feel any buoyancy forces. 
This should be a tolerable assumption during the merger 
itself where the motion of the radio ghosts is dominated by 
advection by the ambient gas. 

In our simulation we store the position, pressure and 
density for all gas particles as a function of time. We as- 
sume that each particle consists, beside the gas, of a certain 
amount of non-thermal radio plasma. Thus each particle rep- 
resents gas whose evolution is calculated by the simulation 
and radio plasma whose evolution is derived from the gas 
properties. For the starting time to before the merger, we 
assign all particles, the same power-law electron momentum 
distribution. The local pressure defines the local magnetic 
field strength due to the fixed ratio Pb/P^os and the com- 
pression ratio, see Eq. Q . We assume that the magnetic field 
within the rarefied plasma bubbles is tangled on small scales 
and can be approximated, together with the relativistic par- 
ticles, by a 7 = 4/3 equation of state. From to the recorded 
pressure history of each particle alone, we can calculate the 
momentum distribution of electrons in the radio plasma, see 
Eq. (|1J . The luminosity of the radio plasma is computed for 
an observing frequency of = 100 MHz and an initial spec- 
tral index of a = 2.5 usin g the standard integration kernel 
for synchrotron radiation (iRvbicki fc Lightmaru.lQTgj) . 



4 RESULTS AND DISCUSSION 

The merger produces shock waves that propagate in both 
directions along the line that connects the centres of the ini- 
tial clusters. While the shock generates only a small jump 
in density, the temperature in the shocked regions is about 
one order of magnitude above that in regions in front of the 
shock (see Fig. . The effects of the merger on the radio 




Figure 3. The projected 'potential' radio luminosities for 
1.13 Gyr old radio plasma, where Ps/Pgas = 0.01. The radio 
plasma is identically to the gas distributed, see Sec. |1] For com- 
parison the bolometric surface X-ray luminosity, Lx = 1-2 X 
10^4ergs-i mgas/(/imp) Pi/(/^rnp) (kTi/keV)i/2 iEke et alj 
\l99ft) . is given. Contours are at , lO*^, 10*3, lO** and 
10'*® ergs~^ h^ Mpc~^. The total bolometric X-ray of the clus- 
ter is 2 X 10'*'*ergs~* and the emission-weighted temperature is 
3keV. 

plasma can be seen in Fig. |21 which shows the radio lumi- 
nosity in a slice through the cluster. The slice was made 
along the plane which contains the both centres-of-mass of 
the initial two clusters. The most remarkable structure is the 
prominent ring-like feature with a diameter of about 1 Mpc 
(see Fig. I^J 1.13 Gyr after the release of the radio plasma. 
This structure corresponds to the outgoing shock waves seen 
in the temperature (see Fig.Qlower panel). One can clearly 
see the flaring of radio plasma at the two outgoing merger 
shock waves. It is apparent that after about 1 Gyr the merger 
shock waves can still revive the fossil radio plasma. In con- 
trast, it is striking that the cluster centre is virtually void of 
any luminous sources. Prerequisite for a reanimation of the 
plasma is a sufficiently low magnetic field strength. Only if 
the strength is as low Ps/Pgas ~ 1% revived structures can 
be seen. Higher magnetic field strengths result in too fast 
an ageing of the plasma, such that the shock waves cannot 
revive the plasma. Furthermore, in the early stage of the 
merger, when the shock waves pass through the centre of 
the cluster after about 0.5 Gyr, the even younger plasma 
does not flare up significantly. 

The much higher occurrence of radio relics at peripheral 
locations in the cluster is a result of two factors: In the centre 
of the cluster, i.e. within a region a few hundred kpc in diam- 
eter, the radio plasma ages much faster because the pressure 
is much higher and, thus, also the (assumed) magnetic field. 
This causes considerably higher radiation losses. In the in- 
ner region of the cluster the luminosity of the plasma dies 
down after ~ 0.5 Gyr, whereas in the periphery the luminos- 
ity decreases much more slowly. The second reason is related 
to the shock compression. As the shocks sweep through the 
cluster, their strength varies with the ICM sound speed. The 
shock waves are relatively weak in the hotter cluster centre 
but steepen when they pass the cooler, outer regions of the 
cluster. Therefore, the compression factor of the shock in- 
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creases and, thus, the abihty to revive radio ghosts. Under 
the condition that the ratio Ps/Pgas is as low as « 1%, 
the two shock waves can hghten up fossil radio plasma in 
the outer regions of the cluster which is about 1 Gyr old. 
In stronger fields, synchrotron losses are too severe and in 
weaker fields, higher energy electrons are required to emit 
at the observing frequency, which have more severe inverse- 
Compton losses. For Ps/Pgas ~ 1% the two losses roughly 
match. We should point out again that Fig. Inland |H] are lu- 
minosity probability maps and not actual maps. Therefore, 
we do not expect radio relics over the entire shock surface. 

We find additional bright spots in the periphery of the 
cluster which are not related to the shock fronts. These spots 
exist even if Ps/Pgas is as high as « 10%. Following the evo- 
lution of merger, one can notice that these spots appear in 
regions where gas from the periphery is flowing into the clus- 
ter. Since the radio plasma from regions that lie further out 
has aged less than more centrally located plasma, a merely 
moderate compression may cause a noticeable increase in 
the luminosity. This provides an additional mechanism to 
generate radio relics, which again only takes place in clus- 
ters that are not in equilibrium. Both mechanisms, gas infall 
and merger shock compression, strongly favour radio ghosts 
that are located at least at a distances of a few hundred kpc 
from the cluster centre. 

We now compute the radio surface brightness. For this 
purpose we convolve the luminosity of radio plasma with a 
hypothetical distribution of radio plasma. Here, we assume 
that the distribution of plasma follows the distribution of 
gas. This approach disregards the fact that there are - prob- 
ably - only few individual plasma object in a cluster, but 
we interpret the distribution of the entire cluster as a prob- 
ability that plasma is located at a certain position. We can 
now infer the radio luminosity of the cluster. Given that the 
distribution of plasma reflects the spatial probability distri- 
bution of radio ghosts, the luminosity corresponds to the 
probability to flnd radio relics at a certain position. We find 
that by far most of the luminosity projected onto a plane 
perpendicular to the shock fronts comes from the shock re- 
gions (see Fig. |^. Even the projection leaves the central 
region of the cluster virtually free of emission. This result 
indicates, that if the distribution of radio ghosts follows the 
distribution of the gas, i.e. if ghosts are most likely found in 
the centre, relics are expected to be observed almost exclu- 
sively at the location of the shock fronts. 

Upcoming radio telescopes, such as GMRT, LOFAR 
and ALMA, will identify more radio relics and will measure 
their spatial distribution. They will be able to verify whether 
relics are predominantly found along outgoing merger shock 
waves at di stances of a few hundred kpc from the centre of 
the cluster jEni31in fc Briiggenll2002l) . Thus, in the near fu- 
ture, observations of a representative sample of radio relics 
will put our model to test. 

In summary, on the basis of cosmological high- 
resolution simulations of a galaxy cluster merger with two 
progenitors of the mass ~ 1.6 x 10^"^ Mq we have 
been able to calculate, both, the fading of non-thermal 
radio plasma and its re-energisation by shock waves. It was 
found that cluster-wide shock fronts are capable of reviving 
~ 1 Gyr old radio ghosts if the ratio Ps/Pgas is as low as 
1%. Magnetic field strengths of this order are not unlikely 
in the cluster periphery. Finally, the vast majority of radio 



relics are expected to be located at typical distances of a 
few hundred kpc from the cluster centre. 
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